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Efficiency of electron heating through microinstabilities generated in the transition region of a 
quasi-perpendicular shock for wide range of Mach numbers is investigated by utilizing PIC (Particle- 
In-Cell) simulation and model analyses. In the model analyses saturation levels of effective electron 
temperature as a result of microinstabilities are estimated from an extended quasilinear (trapping) 
analysis for relatively low (high) Mach number shocks. Here, MTSI (modified two-stream instability) 
is assumed to become dominant in low Mach number regime, while BI (Buneman instability) to 
become dominant in high Mach number regime, respectively. It is revealed that Mach number 
dependence of the effective electron temperature in the MTSI dominant case is essentially different 
from that in the BI dominant case. The effective electron temperature through the MTSI does not 
depend much on the Mach number, although that through the BI increases with the Mach number 
as in the past studies. The results are confirmed to be consistent with the PIC simulations both in 
qualitative and quantitative levels. The model analyses predict that a critical Mach number above 
which steep rise of electron heating rate occurs may arise at the Mach number of a few tens. 

PACS numbers: 52.50.Lp, 52.35.-g 



I. INTRODUCTION 

Efficiency of electron heating in high Mach number col- 
Usionless shocks is one of the outstanding issues of space 
plasma physics as well as astrophysics. It is known from 
observations that quite efficient heating occurs in ex- 
tremely high Mach number shocks like SNR (super nova 
remnant) shocksi"— , while in-situ observations of earth's 
bow shock mostly show inefficient heating of electrons^. 
This is probably due to the fact that a dominant elec- 
tron heating process and its efficiency through a variety 
of microinstabilities in shock transition regions strongly 
depend on Mach numbers. Typical Mach number of 
the earth's bow shock is < 10, while that of the SNR 
shocks is of the order of 100 - 1000. In the late 1980s 
a new process of strong and rapid electron heating in 
extremely high Mach number shocks was proposed^ii. 
Refs. 6 and 7 discussed the so-called two step insta- 
bilities generated by a high velocity reflected ion beam 
which destabilizes BI (Buneman instability) followed by 
ion acoustic instability. A necessary condition for the 
BI to be destabilized in the foot of quasi-perpendicular 
shocks is MA/yj% > 20 ~ 3(^^, where Ma is the Alfven 
Mach number and (3^ = SnnTe/B^ the electron plasma 
beta. If the above condition is satisfied, rapid electron 
heating with saturation temperature cx occurs and 
the ion acoustic instability sets in leading to further elec- 
tron heating^i. This condition is easily satisfied in SNR 
shocks, although it is usually not in the earth's bow 
shock. 

The two step instabilities were again shed light on after 
Ref. 9 by utilizing electromagnetic full PIC (particle-in- 
cell) simulations. They improved understandings of the 
process not only in terms of electron heating but also 
in terms of a production process of non-thermal elec- 



trons, and promoted a number of subsequent simulation 
studiesiSr^. As a result, understandings of electron heat- 
ing or acceleration processes initiated by the BI in the 
context of the shock physics have been extensively devel- 
oped. 

On the other hand, it has been well-known that var- 
ious microinstabilities are possible to get excited also 
in relatively low Mach number supercritical shocks^^^^i. 
Since most of the microinstabilities are inseparable with 
electron dynamics, their nonlinear evolutions in self- 
consistently reproduced shock structures in PIC simula- 
tions are studied only recently^^"— except for Refs. 29- 
31. It is revealed that MTSI (modified two-stream insta- 
bility) becomes dominant in the foot of the earth's bow 
shock or interplanetary shocks at ^ 1 AU^"—. Sim- 
ulation studies on nonlinear evolutions of the MTSI in 
periodic systems were extensively studied in the 1970s 
and 1980s^"— , although these simulations imposed some 
strong assumptions that the code was electrostatio '^^i'^^ 
and unmagnetized ions with small ion-to-electron mass 
ratio were used'^^"'^^. Recently, some developement with- 
out these assumptions are made^ii^^. Ref. 8 showed in a 
one-dimensional simulation with realistic mass ratio that 
long time evolution of the MTSI results in lower cas- 
cade of wave spectrum and associated electron heating. 
Furthermore, it is indicated in two-dimensional simula- 
tion that the MTSI can finally survive even if other pos- 
sible instability like electron cyclotron drift instability 
is present^. All these simulation studies show electron 
heating or acceleration in the long time evolution of the 
system. However, as mentioned already, electron heating 
in the earth's bow shock is seldom observed in-situ. 

Such a discripancy may arise because of lack of sys- 
tematic estimate of electron heating rate in simulation 
studies with more realistic parameters. For example, all 
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of the above simulation studies on the MTSI assume very 
low values of squared ratio of electron plasma to cyclotron 
frequencies, r = tUp^/ill < 10. Importance of r in asso- 
ciation with electron heating in the long time evolution 
of the MTSI is unresolved, while r is known to be a 
crucial parameter controlling nonlinear electron heating 
and acceleration processes in the BI dominant systemic. 
Electron plasma beta is also rather small in most of the 
previous studies (/3e < 0.1). This may exaggerate ef- 
fects of electron trapping in the nonlinear stage of the 
]y[rpgj8|25-_27|33^ The electron plasma beta in the solar 
wind is usually several times higher so that importance 
of the trapping process may be blurred. Hence, in realis- 
tic situations expected to be achieved in the typical solar 
wind conditions an alternative electron heating mecha- 
nism and its efficiency are ought to be considered. 



II. EXTENDED QUASILINEAR ANALYSIS ON 
ELECTRON HEATING THROUGH MTSI 

A. System Configurations 

A foot region of a quasiperpendicular shock is mod- 
eled by a plasma composed of incoming electrons at rest, 
incoming ions, and specularly reflected ions with local 
approximation. The ambient magnetic field is along the 
z-axis, while the incoming and the reflected ion beams 
are streaming parallel and antiparallel to a wavevector 
which is defined in the x — z plane and parallel to the 
shock normal. This is essentially the same system dis- 
cussed in the linear analysis of Ref. 8. Since the MTSI 
grows much faster than ion cyclotron period, ions are 
assumed to be unmagnetized. 



B. Governing Equations 



In this paper, first, electron heating through the 
MTSI in transition regions of high Mach number quasi- 
perpendicular shocks is reconsidered in section II. Here, 
quasilinear diffusion is assumed as an alternative elec- 
tron heating process. The analysis assumes a broad wave 
spectrum and small wave amplitudes which are likely 
to be satisfied in the past simulation studies mentioned 
above2>i23,. An additional assumption that relative phases 
among different wave modes are random is also imposed, 
while this point is currently not evident. Under these as- 
sumptions, taking second order velocity moments of the 
so-called quasilinear equations provides evolution equa- 
tions of electron kinetic energies. These equations are 
numerically integrated by imposing the quasilinear as- 
sumptions in every instantaneous time steps to obtain 
saturation kinetic energies as a fuction of the Mach num- 
ber of shocks. This approach allows to discuss long time 
evolution of macroscopic quantities including the elec- 
tron temperature without any restrictions for parame- 
ters. It is intriguing whether the Mach number depen- 
dence of electron heating throug the MTSI in relatively 
low Mach number regime is different from that through 
the BI in extremely high Mach number regime. If that 
is the case, one may expect presence of a critical Mach 
number above which electron heating rate runs up due 
to the switching of the dominant microinstability^^. In 
section HI, the results are compared with ID PIC sim- 
ulations in which the MTSI is dominantly generated in 
the foot and effective electron temperature just behind of 
the shock is measured. Furthermore, post-shock effective 
electron temperature is measured also in cases that the 
BI gets excited in the foot. In those cases a rough esti- 
mate of the electron temperature is given by using the 
well-known trapping theory. Finally, a summary and dis- 
cussions including estimate of the critical Mach number 
are given in section IV. 



Saturation levels of electron kinetic energy through the 
MTSI are examined as follows. In a Vlasov equation, 

variables are expanded around zero-th order slowly vary- 
ing quantities (with subscript 0) as 



fj - Pjo + fji, 
B Bo + Bi, 
E = El, 



(2) 



where /j, qj, and rrij are a distribution function, elec- 
tric charge, and mass of species j, c is speed of light, 
B and E represent magnetic and electric fields, respec- 
tively. The variables with subscript 1 denote fluctuations 
corresponding to linear oscillations. With a standard ap- 
proach of quasilinear analysis, an evolution equation of 
FjQ is obtained as foUows^i^. 



— — = hm — - y 

n— — CO 



ii;_L^(w.r,k - k\\v\\ - n£t)\^n^^^v^LFQ (3) 

Here, V is the volume of integration, i7 = qB^jmc the 
cyclotron frequency, f || and uj^ the velocity components 
parallel and perpendicular to Bq, k the wavevector which 
is also composed of parallel and perpendicular compo- 
nents (fc||,fcj^), and w,-,k indicates real wave frequency 
which is a function of k and satisfies the linear dis- 
persion relation. The subscript j has been eliminated 
for simplicity. Only resonant wave-particle interactions 
corresponding to zero arguments of the delta function, 
i^r.k — fcjjt'ii — nSl = (rt = 0, ±1, ±2, • • • ), are taken into 
account. Moreover, 
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(5) 



where ^[| = Eiz,i^/Ei^,k, A± = iEiy,i^/ Ei^^]^, Efy^ = 
Eix,k + Ejy^w + ^L,k' J„(C) represents the n-th or- 
der Bessel function with argument C — k±v±/fl. When 
the above expression of ^'n,k is derived, the hnear disper- 
sion relation is taken into account. Note that in eqs.Q- 
([5]) parallel and perpendicular directions are defined in 
terms of Bq. Although this is true for electrons, fur- 
ther note is needed for ions. Because ions are assumed 
to be unmagnetized, distinction of parallel and perpen- 
dicular directions in this sense is meaningless. However, 
their evolution equations can be formally written in a 
similar form by defining parallel and perpendicular direc- 
tions in terms of the ion beam velocity. This corresponds 
to a transformation for the rotation of the susceptibility 
tensor—. 

To discuss the extended quasilinear evolution of kinetic 
energies through the MTSI, we examine parallel and per- 
pendicular kinetic energy densities defined by 
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(6) 



(7) 



After some algebra, the following evolution equations of 
the kinetic energies are obtained. 



dKn 



dt 




UJrM\k\\ I 



(8) 




(9) 



Here, v = w^.k ~ ^n'^o ~ i^^, and ()5=o means a corre- 
sponding value at which the argument of the delta func- 
tion is zero. The spectral energy density of the electric 
field fiuctuations evolves as 



dEl 
dt 



27k£^, 



kJ 



(10) 



where 7k (i) denotes the linear growth rate of the wave 
mode having wavevector k. Conservation of total energy 



is written as 



— ^.0 



K 



V ^ = const. (11) 

k 



Here, rij and Vjo(t) denote density and bulk velocity of 
species j. Note that dK±/dt — for unmagnetized ions 
and Veo = 0. Further, ujr,k and 7k are defined at an 
instantaneous time. Eqs.([5|)- pT|) form a closed set of 
extended quasilinear evolution equations of the system. 
We solve the above set of equations numerically in the 
following subsection. 

Similar analyses were performed by many 
authorsHiH^iS. Ref. 38 introduced detailed proce- 
dure of the method and applied it to estimate of 
saturation of electromagnetic ion cyclotron instability 
driven by ion temperature anisotropy. Refs. 39 and 
40 directly solved the quasilinear equation to discuss 
electron heating and formation of a high energy ion 
tail through the ion acoustic instability. Ref. 33, on 
the other hand, discussed a problem similar to here 
on the MTSI and gave detailed comparisons of their 
results with Id and 2d PIC simulations. However, 
their analysis includes only nonresonant wave-particle 
interactions which in turn only slight electron heating 
which is much less effective compared with a subsequent 
heating process through nonlinear trapping observed 
in their PIC simulation. Furthermore, they imposed 
an electrostatic assumption on whole their analyses 
as well as on their PIC simulations. As they pointed 
out, the electrostatic approximation is valid only when 
Vjo/i^A < (1 + /3e)^^^, where va denotes the Alfven 
velocity, and /?e is the electron plasma beta (ratio of 
electron thermal to magnetic pressures). Usually in 
solar wind f3e is of the order of 0.1 — 1 so that the above 
condition reads vjq/va ^ 1 which is probably satisfied 
only in transition regions of rather low Mach number 
shocks. Since our attention is paid for wide range of 
Mach numbers, we exclude such an assumption. 

Before solving eqs.([8|)- (fTT|) numerically, a few assump- 
tions are imposed. The distribution functions of each 
species remain in (shifted) Maxwellian all the time after 
Ref. 38. By this assumption, it is easy to obtain ujr,k.it) 
and 7k (i) at each time step as instantaneous solutions 
of the linear dispersion relation. The actual distribution 
functions are given by 



4. 



exp 



(12) 

^/Tj\\/mj and vtj± = 



(2^)3/2«,,.,[«; 

for species j. Here Vtj\ 
y^Tj±/mj denote parallel and perpendicular thermal ve- 
locities. As mentioned before, note that the parallel and 
perpendicular directions are difined in association with 
Bq for electrons, while with vjo for incoming and re- 
fiected ions. The bulk velocities of the incoming and 
refiected ions (subscripts j = i and r, respectively) sat- 
isfy zero current condition so that UiVio + rirVro = 0. It 
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is also assumed that the response of the ions to the wave 
fields is electrostatic, while the wave-electron interactions 
are treated as fully electromagnetic. Therefore, the ion 
heating is thought to be parallel to Vio or Vbo- 

Under these assumptions, we obtain the following evo- 
lution equations of kinetic energies. 



dK, 



dt 



lim > 

n 



dt 



lim > 

n 




UJrM\H\ 

dvx_v\{\^n?FeQ)s^O 



(13) 




UJrM\k\\\ 



(14) 



dt 



lim 

y-i-oo TO," 



dv±vl{\Mz^O)\^FM)s=o 



(15) 



Here, Ve — Wr,k + nfle and i>i — cjr,k — kviQ, respectively. 
In this analysis only the linear resonant wave-praticle in- 
teractions are taken into account. Therefore, change of 
kinetic energy of the refiected ions {K.^^^ ) can be neglected 
when we consider the MTSI based on electron-incoming 
ion interactions. That is, here, the refiected ions are as- 
sumed to be a background component satisfying current 
and charge neutralities. Reactions of the refiected ions 
may be essential when influences of the MTSI on the ref- 
ormation process are consideredSi,. However, it is out of 
scope in this paper. In the following we solve eqs. (ll3p - 
(113, (Uni), and 



Tiinii 



K,n + K, 



3, 

k 



const. (16) 



instead of eq. pTj) . Eq. (fT6| is actually used to determine 
ufo- Note that contributions from the terms |n| > 4 for 
J2n eqs. ([T^ and ^(T^ are neglected. The validity of 
such an assumption will be checked in the following sub- 
section. 



C. Numerical Solutions 

First, let us overview in Fig[T] time evolutions of the 
variables calculated in the system for a typical case. 
The initial parameters are /i = mi/m^ = 1836, r = 
uljnl = 104, /3 = 8^ne(Te + Ti)/Bl = + = 




Hit 



FIG. 1. Time evolution of the system. 



0.4, a = rir/n^ = l/3,6'Bfc(= Qbu) = 85.5°, and Ma = 
{vio/vA){^ + Oi)/2a = 6. Here, the initial temperatures 
are isotropic and Pe = Pi- Ma and a determine the 
initial Wio from the zero current condition. 0BkiQBn) de- 
notes the angle between the ambient magnetic field, Bq, 
and the wavevector (the shock normal), k(n). The above 
parameters are fixed throughout the analysis if not spec- 
ified. The field energy is given as appropriately small 
value: E^/Bq = 10~^'^. For the moment, kc/ojpe is 
fixed to unity just for simplicity. The field energy ex- 
ponentially grows in the beginning of the run, although 
other variables have not been much changed in a vis- 
ible manner until the field energy reaches at a certain 
level. Then Vio gradually decreases, while ife|| and Ki 
increase {^lit > 0.2). At the same time, the wave growth 
slows down and it starts damping after fi^t ~ 0.45. Cor- 
respondingly, changing rates of other variables decrease 
and those variables approach some constant values when 
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the field energy becomes sufficiently small. Such a state is 
regarded as a saturation state and those constant values 
are defined as saturation levels. The saturation state in 
this system is achieved when the field energy is damped 
out. In the analysis we check whether the system sat- 
urates basically by eye. However, in some special cases 
where the kinetic energies have not reached at constant 
values for long time the saturation levels are defined as 
the values at t — to + 2il^^ where to is the time at which 
ViQ decreases by 1% from its initial value. Since 2f2~^ is 
thought to be a typical time scale of shock reformation, 
considering the system evolution longer than this time 
scale may not make sense. 

In FiglT]iirej_ keeps almost constant incontrast to ifgH 
or Ki. This reflects which kinetic effects dominantly 
work. The linear growth rate, 7^, can be written by 
linear superposition of the terms such as 



?joexp(-^|„), 



(17) 



where ^en = {^^r + nfle)/ki\Vte\\ and^io = {^r-kvio)/kvu, 
respectively. When ~ 1, contribution from the 

corresponding term to the growth rate becomes non- 
negligible. In this sense S,eOi ^ei, and can be indicators 
how strong the effects of Landau and cyclotron damp- 
ings of electrons, and inverse ion Landau damping are, 
respectively. The bottom panel of FigH] represents val- 
ues of these indicators. It is now clear that the electron 
Landau damping is the strongest kinetic effect and the 
electron cyclotron damping hardly works. The inverse 
ion Landau damping can also work, but may not be so 
strong as the electron Landau damping for this particu- 
lar case. This can explain the fact that increase of Kf.\\ 
is most remarkable, while that of Ke± is unrecognized. 
Contributions from other electron cyclotron harmonics 
{£.en with n = —1, ±2, ±3, ±4) are further ineffective, and 
this is confirmed to be true for wide range of parameters 
(not shown). Therefore, effects of electron cyclotron in- 
teractions will be neglected in the remaining of the anal- 
ysis. The assumption that kc/tUpe = 1 is fixed is removed 
hereafter. 

The upper panel of FigH] denotes saturation lev- 
els of normalized electron temperature, /3e||, sat = 
8TrneTf,\\ sat/ Bq as a function of the initial /Sq — /3oe + f^oi 
when r = 10^. In calculating the right hand sides of 
eas. p^ - ([T5|) . we chose interval of integration with re- 
spect to fc as < kc/ujpe < 3 and grid size in k as 
Akc/ujpe = 0.06. 9Bk = 85.5° is again fixed here. In 
both upper and lower panels the solid line with circles, 
the dashed line with squares, and the dotted line with 
triangles correspond to the cases for /3oe||//3oi = 1,3, and 
1/3, respectively. I3e\\.sat is essentially independent from 
/3o for l3o > 0.4 indicating that heating process is similar 
for relatively high /3o- In the low /3o, on the other hand, 
there is a trend that electron heating is suppressed for 
/3oe||//3oi = 1/3. For sufficiently low /3oe||j number of reso- 
nant particles is very little so that it takes extremely long 
time to heat electrons. Then after the very long time, 
heating is suddenly triggered when field energy becomes 
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FIG. 2. Saturation level of normalized parallel electron tem- 
peratures, as a function of (upper panel) initial /?o and 
(lower panel) r. The three lines correspond to cases with 
different initial temperature ratio. The solid line with cir- 
cles denotes the case with /3oe||//Soj = 1, the dashed line 
with squares /3oe||//So! = 3, and the dotted line with trian- 
gles Poe\\/Poi = 1/3, respectively. 



large enough. Such an unnatural time evolution occurs 
because of the assumptions that only the resonant wave- 
particle interactions are taken into account and that the 
distribution functions are always Maxellian. Indeed, it is 
known from PIC simulations that effects of particle trap- 
ping, which have not been included in this analysis, be- 
come important for low /3o case o^'^^ . The initial electron 
to ion temperature ratio, l3oe\\/Poii affects the saturation 
levels. For large (small) /3oe||//3oj(= 3(1/3)), Pelsat is 
high (low). Comparing with the case of Poe\\/l3oi — 1, 
l^iol gets closer to unity because vu gets larger when 
/3oe||//3oi becomes 1/3. Hence, the inverse ion Landau 
damping works more effectively so that electron heat- 
ing becomes less dominant. The r dependence of /3e||_ sat 
is shown in the lower panel of Figl2] for /?o — 0.4. As 
expected from the linear analysis^, dependence on r is 
rather weak. This implies that using small r in simula- 
tions does not lead to unrealisitic results for long time 
evolution of the MTSI. 

Figl3] represents Ma dependense for r = 10^ and 
/3o = 0.4. Three black lines (the solid line with circle, 
the dashed line with squares, and the dotted line with 
triangles) in the upper panel denote the maximum sat- 



6 



10^ 




10 

Ma 

FIG. 3. Mach number dependence of (upper panel) maximum 
saturation electron temperatures, and (lower panel) cor- 
responding wave propagation angles, dsk- The three black 
lines correspond to cases with different initial temperature ra- 
tio. Their definition is given in the caption of Fig[2] The thick 
gray solid (dashed) line is obtained from eg. ljlSp (ea. (|19p V 



uration levels of /3e||,sot in varying Osk again for differ- 
ent initial temperature ratios (/3oe||//3oi — 1,3, 1/3). The 
corresponding values of Obu is indicated in the bottom 
panel. /3e||, sat does not increase so much, even though 
the Ma gets larger. As already mentioned, the system 
saturates when available field energy is lost and no more 
waves grow. Such a state may be achieved when effect 
of electron Landau damping becomes efficient. If it is 
the case, ry^eo ~ 1 should be satisfied at the saturation 
state. Here 77 is a factor of the order of unity and is em- 
pirically determined later. By using Uj^y = Tgn/me and 
Lo/k\\ w Viq/ cos 9 Bk, 

The gray thick solid line in Figl2] is obtained by sub- 
stituting the values of Osk of the solid line in the bot- 
tom panel into ea.([T8|. Here, rj = 2.6 has been as- 
sumed. Although dependence on initial tempetarure ra- 
tio, /3oe||//3oi, is not trivial in eq. ()18p . this may give a rea- 
sonable order of magnitude estimate. Remember that in 
the MTSI a cross-field ion beam interacts with oblique 
whistler waves. The whistler waves which can interact 
with the high speed ion beam should have appropriately 
high phase velocities. The phase velocity of an oblique 
whistler wave is proportional to cosOsk^ so that a high 
Mach number shock generates high phase velocity, i.e., 
less oblique, whistler waves. This effect appears in the 
bottom panel. As a result. Ma/ cos Osk does not change 



TABLE I. Injection Plasma Parameters for Simulation 





MAin 


T 


mi/rrie 




ft 




Run A 


4 


16 


625 


0.3 


0.1 


84 


Run B 


6 


16 


625 


0.3 


0.1 


81 


Run C 


10 


16 


625 


0.3 


0.1 


79 


Run D 


8 


100 


64 


0.3 


0.1 


90 


Run E 


16 


100 


64 


0.3 


0.1 


90 


Run F 


30 


100 


64 


0.3 


0.1 


90 



so much even the Mach number varies, so that /3e||. sat 
remains in the same order for wide Mach number regime 
examined here. Only the dotted line with triangles in the 
upper panel is weakely decreasing function of the Mach 
number. In this relatively high initial ion temperature 
case, as mentioned earlier, ion kinetic effect seems to be 
non-negligible. As a result, ion heating is more remark- 
able than in the cases with ^oe|| /Poi = 1,3. Furthermore, 
because of small growth rate due also to the ion kinetic ef- 
fect, it takes rather long time for the system to approach 
a final saturation state. Therefore, the saturation levels 
were actually estimated at i = to + 2^2"^ as mentioned 
before. 

Note that the condition i^eo ~ 1 gives essentially same 
expression given in Ref. 33 which derives saturation en- 
ergy due to electron trapping. Hence, the expression in 
eq. (fT8|) may be applicable in also the low beta case for 
MTSI. However, because Ref. 33 assumed that ^ cos^ Osk 
is constant throughout their analyses, their saturation 
temperature increases with M\. This might be true if the 
Mach number is sufficiently small and the electrostatic 
approximation is valid. But it is now clear that constant 
/icos^ Osk assumption is invalid at least for Ma > 5. 

III. ELECTRON HEATING REPRODUCED BY 
ID PIC SIMULATION 

Here, shock waves are actually reproduced by using 
one-dimensional PIC code and the resultant electron 
heating through MTSI and BI in the transition region 
is discussed. A shock is produced in the simulation 
domain by using the so-called reflecting wall method. 
An upstream plasma consisting of equal numbers of 
ions and electrons is continuously injected from the left- 
hand boundary (x — 0) and carries a uniform magnetic 
field, Bo — {Box,0, Bqz), and a motional electric field, 
Ey — UinBoz/c, where Uin is the injection flow velocity in 
the positive x-direction. At the right-hand boundary the 
particles are reflected. As a result of the interaction be- 
tween the incoming and reflected particles a shock wave 
is produced and propagates in the negative a;-direction. 
Thus the simulation frame is the downstream rest frame. 
The grid size is set to be Ax « 0.83A_De and 200 par- 
ticles/cell are distributed for both ions and electrons at 
the injection boundary, where Xoe denotes the electron 
Debye length. Physical parameters are summarized in 
TABLE n For runs A, B, and C (D, E, and F), MTSI 
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FIG. 4. Electron heating observed in one dimensional PIC 
simulations for (upper panels) Run A and (lower panels) Run 
F. The left panels show electron phase space distributions 
in (upper left panel) Vez — k and (lower left panel) Vex — 
X. The gray solid lines denote profiles of magnetic field Bz 
component. The right panels indicate electron distribution 
functions in v^.^ and integrated over space corresponding 
to the arrowed regions in the left panels. 
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FIG. 5. Comparison of normalized electron temperatures 
between the PIC simulations and the model analyses. The 
markers labeled A ~ F are estimated as averaged effective 
electron temperatures downstream of the shocks reproduced 
in the PIC simulations. The dashed lines are obtained from 
the model analyses, ea. (|18|) for A ~ C and eQ (|19p for D ~ F, 
respectively. 



(BI) gets excited in the foot during each reformation cy- 
cle. 

Fig|3] shows two examples of snap shots of electron 
phase space — x at fi^i — 6.02 for Run A (upper left 
panel) and Vx — x a,t flit — 8.67 for Run F (lower left 
panel). The gray solid lines denote profiles of magnetic 
field Bz component. In both runs velocity distribution 
functions corresponding to the downstream areas indi- 
cated by arrows in the left panels are represented in the 
right panels. Nonadiabatic electron heating downstream 
of the shock is clearly observed. Parallel heating is su- 
perior to the adiabatic perpendicular heating in Run A, 
while strong perpendicular heating, or heating parallel 
to the beam velocity, is seen in Run F. The effective 
temperatures, K'^-^^ , just downstream of the shocks are 
estimated by using the same definition as eqs.iH]) and 
(O. We averaged over the particles distributed in the 
region of Xos ^ x < Xos + Pi, where Xos is the position of 
the magnetic overshoot and pi the ion cyclotron radius 
defined by upstream flow velocity and ion cyclotron fre- 
quency. Because of nonstationarity of the shock, K'^^^ is 
not constant in time. Therefore, K'^^^ is also averaged 
in time a few shock reformation periods. The results are 
plotted in Figl5]as a function of the injection Mach num- 
ber, MAin- Here, K*^^^ is normalized to the upstream 



magnetic pressure, i?oi/ 



?7r, and rewritten as (3'^^^ . For 



Runs A, B, and C, is almost independent on MAim 
which is consistent with the extended quasilinear analy- 
sis shown in the previous section. The dashed line just 
below the points A to C is drawn by using eq. p^ with 
assuming a = 1/3, rj = \, and Ma = MAin- 

On the other hand, (i^l^ efficiently increases with Maiu 



for Runs D, E, and F, where BI gets excited in the foot. 
Nonlinear development of the BI was well studied in the 
past. Especially, Ref. 10 analyzed it in detail by per- 
forming one dimensional PIC simulations with periodic 
boundary conditions modeling a part of the foot region 
in a perpendicular shock geometry. They showed that 
when r ^ 1, the BI produces well defined electron holes 
in its early nonlinear stage and some additional heating 
or acceleration occurs in the course of further long time 
evolution. Although a saturation level of the instability 
is reduced a little in two dimensional system, electron 
trapping still plays an important role through oblique 
modes^S. It is infered from these results that size of the 
electron holes in phase space may give a lower limit of 
P'j.l^ . The amplitude of an electron hole was estimated by 
assuming that electrostatic wave energy density is com- 
parable with drift energy density of electrons in the ion's 
rest framei^i^. For our case, it reads 



ex 



Ml 



1 



a /I 



7/6' 



(19) 



where we assume the BI based on electron-reflected ion 
interactions which usually grows faster than the BI based 
on electron-incoming ion interactions so that I^u/va — 
2Ma/ (l + ck) is used as the relative drift velocity. As seen 
in the above expression, fi^l^ increases with M\ which is 
proportional to the system free energy. This functional 
dependence is also derived in Refs. 6, 7, and 11. One 
can confirm that eq. p9)) may give a good estimate of the 
lower limit of jS'^l^ by plotting the corresponding values 
in FiglS]for a = 0.5 and = 64 as the dashed line just 
below the points D to F. One should note that Ma in 
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eq. (IT9l) is not necessarily coincide with Maiu in Figl5] 
Because the indicates the local Mach number of the 
foot in the shock rest frame, while the Maiu denotes the 
injection Mach number in the simulation frame, i.e., the 
downstream rest frame. It is also hard to define the local 
Mach number of the foot from the simulation data, since 
the shock is not time stationary. Therefore, the dashed 
line should be regarded only as a guide. 



IV. DISCUSSIONS 

In the previous section it was confirmed that the ex- 
tended quasilinear analysis for the MTSI and the trap- 
ping analysis for the BI give consistent results with the 
PIC simulations. Here, let us compare the results of these 
two analysis. The thick gray dashed line in Figl3] denotes 
/3||^ obtained from eq.^ for = 1836, a = 0.5. This 
line intersects at Ma ^ 24 with the dotted line with tri- 
angles which is (3e\\,sat for the MTSI with l3oe\\/Poi — 1/3- 
The gray thick solid line based on eq. lfTS)) likely to inter- 
sect with the gray thick dashed line at Ma ~ 45. De- 
pending on parameters, an actual intersection may occur 
in the range of 20 < Ma < 50. If this is written as Af^, 
the dominant electron heating process possibly switches 
from the MTSI for Ma < M^ to the BI for Ma > MX- 
For relatively low Mach numbers like in the earth's bow 
shock, /3e||, sat IS more or less constant with 10°""^ be- 
cause the MTSI becomes dominant. If /3e.sat — 5 and 
T = 10^ are assumed, corresponding electron tempera- 
ture is ~ lOOeV which is consistent with a typical tem- 
perature observed downstream of the earth's bow shock. 
However, addiabatic heating due to increases of magnetic 
field and cross shock potential also results in the down- 
stream electron temperature similar to this value. This 
may be the reason why remarkable electron heating was 
seldom observed in near earth shocks in the past. But 
one may capture the nonadiabatic parallel electron heat- 
ing if upstream electron beta is low enough (/3eo < 0.1), 
since eq. (fT8|) is still valid for rather low beta cases where 
the trapping effects become essential as mentioned in sec- 
tion II. Furthermore, as shown in the upper right panel 
of FigHl local temperature anisotropy, Te\\/Te± > 1, in 
a transition region may be observed as a result of strong 
MTSI in some parameter regimes. On the other hand, 
the effective electron temperature significantly increases 
with being proportinal to in Ma > M\. Although 
this is consistent with the past studies^i^iii, the satura- 
tion electron temperature estimated in this study seems 
to be rather smaller. This should be because of that the 
heating process in highly nonlinear stage including the 
second step ion acoustic instability is neglected. Elec- 
tron temperature seen in Ref. 10 seems to be about one 
order higher than the estimate given by eq. (IT51) . This 
implies that an actual M\ may appear at a little smaller 
value. 

Consistency of the results of the extended quasilinear 
analysis and the PIC simulations for the MTSI dominant 



cases indicates that the resonant wave-particle interac- 
tions are essential in electron heating through the MTSI. 
In other wards, nonresonant wave-particle interactions 
which have been neglected in the analysis may not be 
effective. It should also be noted that using unrealisti- 
cally small r in PIC simulations is justified as far as the 
MTSI is concerned because of the weak dependence on r 
(bottom panel of Fig. 2) in contrast to the BI dominant 
systems. 

The transition at Ma — M\ might not be so drastic 
in a realistic case. In the analyses presented here, all 
other possible candidates of microinstabilities have been 
neglected. For example, ECDI (electron cyclotron drift 
instability) is one candidate-S and it might become im- 
portant around this critical Mach number, although there 
are some negative indications for the ECDI to become 
dominant. For instance, it is known that the saturation 
level of the ECDI is usually not so high--. In addition, 
if the ECDI gets excited simultaneously with the MTSI, 
the MTSI becomes dominant for wide parameter range^^. 
Furthermore, it is shown that the ECDI can be impor- 
tant only when r is of the order of unity Hence, the 
Mach number regime where the ECDI is dominant seems 
not to be so wide, but we do not remove the possibility 
that the ECDI gives some contributions to electron heat- 
ing around M\. Because the ECDI is insensitive to the 
electron Landau damping which may strongly affect to 
the BI in this Mach number regime. Further, it is known 
that the MTSI cannot get excited in extremely high Mach 
number shocks like Ma > 30 ~ 40^.. There are also other 
possible microinstabilities in higher dimensional cases^i. 
Contributions from them should be carefully estimated. 
Some of them have been discussed recently by perform- 
ing two dimensional PIC simulationsi2iiSi2i. The multi- 
dimansionality may give another important contributions 
to electron heating. Ref. 42 pointed out that electrons 
are accelerated in the so-called rippled structure which 
is the ion scale structure along the shock surface. Never- 
theless, what should be emphasized here is that the Mach 
number dependence of the effective electron temperature 
is systematically different between high and low Mach 
number regimes. 

In the present extended quasilinear analysis the damp- 
ing of the wave energy in the late stage (e.g., ^it > 0.45 
in FiglH) may be an artifact. Because of the assump- 
tion that the distribution function is always Maxwellian, 
the so-called plateau of the distribution function is never 
produced in this system. The field energy might saturate 
earlier at a certain level if the plateau is produced. In 
this regard, it might be better that the saturation lev- 
els are defined as the values at the time when the field 
energy becomes the maximum. However, the resultant 
saturation levels based on the two definitions, i.e., the 
saturation temperatures estimated at times correspond- 
ing to the maximum and the later sufficiently small field 
energies, are not so much different from each other. Al- 
though it is possible to solve the quasilinear equation, 
eq.Q, directly as done in Refs. 39 and 40, that is the 
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future work. 

The explorations of the inner heliosphere will be under- 
way through the BepiColombo mission. The perihelion 
point of the mercury is about 0.3AU where quite high ve- 
locity IPSs which have not been decerelated are expected 
to be observed. Rcf. 43 estimated the propagation speed 
of one of the IPSs observed in August 1972 and concluded 



that it approached to ~ 2000 km/s at 0.3AU from the 
sun. More optimistic estimate for the same IPS is given 
in Ref. 44 as 4000 km/s at 0.3AU. These results imply 
that the Mach number of this IPS approached to several 
tens. Therefore, the transition of the electron heating 
efficiency at M^, if present, will be possibly observed in 
sititi^. 
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